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文章 编号 :1005-3085(2010)05-0859-06 


一 类 捕食 模型 椭圆 方程 解 的 稳定 性 * 


查 淑 玲 :2， FH, WMR 
(1- 陕西 师范 大 学 数学 与 信息 科学 学 院 ， 西 安 710062; 
2- 渭南 师范 学 院 数 学 与 信息 科学 系 ， 渭 南 714000; 3- 陕西 科技 大 学 理学 院 ， 西 安 710021) 
摘 E: 本 文 在 Dirichlet 边界 条 件 下 研究 了 一 类 带 种 内 相 食 的 捕食 模型 的 平衡 态 问题 。 利 用 谱 分 析 和 分 
歧 理论 的 方法 ， 给 出 了 发 自 半 平凡 解 的 局 部 分 歧 ， 得 到 了 局 部 分 歧 解 的 结构 。 同 时 ， 利 用 线性 算 
子 的 扰动 理论 ， 证 明了 该 局 部 分 歧 正 解 是 无 条 件 稳 定 的 。 最 后 利用 整体 分 歧 理 论 ， 将 局 部 分 歧 延 
拓 为 整体 分 歧 ， 给 出 了 分 歧 曲 线 随 参 数 的 整体 走向 。 所 得 结果 表明 ， 当 参数 满足 一 定 条 件 时 ， 捕 
食 者 和 被 捕食 者 能 够 在 给 定 区 域内 共存 。 
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分 类 号 : AMS(2000) 35K57 中 图 分 类 号 : 0175.26 文献 标识 码 : A 
1 引言 


生物 界 中 动物 的 种 内 相 食 现象 屡见不鲜 ， 种 内 相 食 是 发 生 在 同一 物种 间 以 本 物种 个 体 为 食 
物 的 现象 ， 对 于 带 种 内 相 食 的 捕食 系统 的 研究 在 文献 中 出 现 较 少 ， 因 此 本 文 研究 文献 [1 带 种 
内 相 食 的 反应 扩散 方程 


= D,Autru(l— k) ow 4) 





ea bnuv nv? 
Ut = D,&v + l+nhu cu l+nhu? 


Sip ARE Laplace HF, u(x,t) Mo(2,t) URR RAMARAMA BE, EBROIN E 
常数 ， 其 生物 意义 参见 文献 [1]。 为 了 减少 系统 (1) 所 含 参数 的 个 数 ， 假 定 扩散 系数 相等 ， 令 


@=nhu, t=(na/r)v, t=(r/a)t, d=ac/r, m=nhK, k=ab/(hr), 


无 量 纲 化 后 的 方程 简化 为 
= Au +au(l - #) - i 
-Av+ BY dv- ah. 


TEE E AE on (2) 图 灵 不 稳定 的 条 件 ， 并 使 用 数值 
模拟 得 到 了 不 同 条 件 下 时 空 动力 系统 的 图 灵 模 式 。 本 文 将 考虑 系统 (2) 的 平衡 态 问题 非 负 解 的 
存在 性 和 稳定 性 。 其 平衡 态 方程 为 


Au- au(l — #)+ $% =0, reg, 


(2) 





-Av - $2 +dv+ $ =0, zeg, (3) 


u=v=O0, z EOR. 
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2 ”准备 知识 


首先 给 出 一 些 相关 的 基本 结论 。 
设 g(z) e C(Q)， 线 性 特征 值 问 题 





—Au+q(z)u=dAu, TEQ, u=0, TEN 


的 所 有 特征 值 可 排列 为 Mt(9) < Ala) < As(q) < …， 其 中 入 (gq) 是 单 重 的 主 特征 值 ， 且 关 
F q(x) 严格 单调 递增 站。 简 记 入 = 和 1(0)。 
考虑 非 线性 边 值 问题 


—-Au=au(l—u/m), «EQ, u=0, rEan. (4) 


由 文献 [3,4] 3, “a> IN, ARE (4) 存在 惟一 正解 ， 不 妨 记 为 6。 易 得 0 < 0 < m。 
利用 极 大 值 原 理 给 出 (3) 非 负 解 的 一 个 先 验 估计 。 
引 理 1 设 (wv) 是 (3) 的 一 个 非 负 解 ， 则 a > Xi HO<u<0, O<v<k(1+m). 
为 方便 叙述 ， 下 面 引入 一 些 记号 : 


C29) = {ue C): ula =0), X= CRM) xA), Y=0°0). 


定义 非 线性 算 子 下 : R x XY, HF 





2 
F(d;u,v) = (- Au-au(1 - =) + er 一 Av 一 了 十 dv 十 Tru): 

因 (3) 在 平凡 解 (0,0) 处 的 线性 化 算 子 的 所 有 特征 值 为 {和 i 一 a}U {A +d} A, WER 

引 理 2 当 a > Xi 时 ，(0,0) 是 不 稳定 的 ， 当 a < Xi 时 (0,0) 是 渐 近 稳定 的 。 

Ha> Xi 时，(3) 有 半 平 凡 解 (6,0)。 设 (3) 在 (6,0) 处 的 线性 化 算 子 为 L， 经 计算 易 得 以 下 
引 理 。 

引 理 3 当 和 Xi(-kb/(1+b)+d > 0 时 ， 半 平凡 解 (9,0) 是 渐 近 稳定 的 ; 当 A (—k0/(1+0))+ 
d < 0 时 ， 半 平凡 解 (9,0) 是 不 稳定 的 。 


3 ”分歧 解 的 存在 性 与 稳定 性 


本 节 以 d 为 分 歧 参 数 ， 讨 论 (3) 发 自 半 平 凡 解 (6, 0) 的 分 歧 。 易 知 (u,v) 是 算 子 方程 F(d; u,v) 
二 0 的 解 当 且 仅 当 (w,v) 是 (3) 的 解 。 

定理 1 设 a > à HA(—kO/(1 +0) < 0， 则 存在 充分 小 的 6 > 0， 使 得 对 任意 d E 
(do 一 ô, do 十 6)，(3) 存在 惟一 的 正解 (w(s),v(s))， 且 满足 0 < u(s) < 9, 0 < v(s) (0< s< ô), 
其 中 do = —A1(—k0/(1 + 8)» 

证 明 设 (3) 在 (d;w,v) = (do;9,0) 处 的 线性 化 算 子 为 Lo， 设 N(Lo) 为 核 空间 ，R(Lo) 为 值 
域 空间 。 因 此 





N(Lo) = span{(¢i,%1)7}, dimN(Lo)=1, 
其 中 (p1, p1) 满足 
| -Ap -ap+ yp = -IY 加 


—Ay + doy — $y =0. 
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p 是 算 子 (-A — k6/(L+0)) 的 主 特征 值 对 应 的 主 特征 函数 ， 且 满足 
widz =1, ol = (—A-a+2a6/m)~*(-641/(1 + 6)). 
Q 


由 广义 极 大 值 原理 知 pi < 0。 
Xt Lo 的 伴随 算 子 [有 NN(L) = span{(0,~1)7}. Abt 


R(Lo) = {C0)7 EX: fona = o}, codimR(Lo) = 1. 


i Li = Fidu w) (doju, v) Ale Li - (p1, Y1) = (0, 41)" Z R(Lo)- 
根据 文献 [四 中 Crandall-Rabinowitz 分 歧 定理 ， 存 在 充分 小 6 > 0 及 01 类 函数 w(s) : (一 6, 8) 
—> R, (p(s), p(s)) : (一 6,6) 一 了， 使 得 当 0 < s < 6 时 
d(s) = do + w(s), uls) =0 + spit+ sp(s), v(s) = spı(s) + sy ls) 
Æ F(d; u,v) = 0 FE (do; 6,0) 附近 的 惟一 正解 ， 即 (uls), v(s)) 是 (3) 的 正解 ， 且 满足 
w(0)=0, yp(0)=0, ¥(0)=0, (p(s),(s)) € R(Lo). 
下 面 给 出 分 歧 解 的 稳定 性 。 由 定理 1 的 证 明 可 知 ，0 是 Lo ie, 
由 文献 [6,7] 知 ， 存 在 两 个 分 别 定 义 在 do 和 0 邻 域内 的 函数 4 一 (y(d), S(d)) M s 一 (n(s), 
T(s)), 18 
(y(do), S(do)) = (0, (¢1,41)7) = (n(0), T(0)), 
并 且 
Fu,v)(d; 0,0)S(d) = y(d)S(d), |d- do| & 1, 
Feu v) (d(s); u(s), v(s)) T(s) = n(s)T(s), O<s«1. 
这 里 
S(d) = (wr (d),we(d))", — T(s) = (Xi(s),Xa(s) ， 
mH. y (do) #0, XE n(s) A 0(|s| <1), W 
sd'(s)7'(do) 
lim, n(s) 
4} IC A (u(s), v(s)) 的 稳定 性 由 w(s) 确 定 。 当 nn(s) > 0 时 ，(w(s),v(s)) 是 稳定 的 ， 当 n(s) < 
0 时 ，(w(s),v(s)) 是 不 稳定 的 。 
首先 确定 (do) 的 符号 。 由 引 理 3 知 y(d) = A1(—k0/(1+0)) +d, iy (do) = 1. 
其 次 确定 d'(s) 的 符号 。 取 d(s) = do + w(s)， 只 需 确定 w'(s) 的 符号 。 将 df(s) 及 分 歧 
解 u(s),v(s) 代入 (3) 的 第 二 个 方程 得 
SeS EN Pt = (6) 
两 端 同 除 以 s， 再 关于 s 在 s = 0 处 微分 ， 两 边 同 乘 以 妨 然后 在 Q 上 积分 ， 并 利用 Green 公 
式 得 


三 一 ]. 





vi (kei — (1+ 0)w1) 
o= 人 (T+ 0 一 一 一 一 一 一 dz. 
由 于 pl <0, yı >0, Mw’(0) < 0，d'(0) < 0。 从 而 当 0 < s。< 6 时 ，d(s) <0, n(s)>0. 

定理 2 设 定理 1 的 条 件 成 立 ， 则 由 定理 1 确定 的 分 歧 解 (u(s),v(s)) (0 < s < 6) 是 渐 近 稳 
定 的 。 


862 I E 数 学 学 R 第 27 郑 


4 全 局 分 歧 
AT MHS AaB, Sw=0-u, x=v, 则 将 0<w<6, 0<X 代 入 (3) 得 


—Aw = aw 一 28w + 795% + Fi (w, x), 


-Ax = -dx + $$x + Fal, x), (7) 


w=x=0. 


记 Fw,x) = (Fi, 把)， 显 然 下 是 连续 函数 ，F(0,0) = 0A Fit Frechét 导数 Fiu (0,0) = 0. 
E Dirichlet 边界 条 件 下 ， 取 KK 为 算 子 (一 人 ) 的 逆 算 子 ， 定 义 


2a 0 
T(d;w,x) = (akw = 元 五 (gw) 十 Ky 


x) + KFiw,x), 
-dKx+K(- Os x) +KF,(w,x)). 


W G(d;w, x) = (wxX) 一 T(d;w,X)， 则 当 0 < w < 0, 0 < x, Gldjw,x) = 0 的 解 就 是 (7) 的 
非 负 解 ， 易 见 G 是 C1 函数 且 G(qd;0,0) = 0， 所 以 w= 0 一 w, v= xÆ (3) HAE. + 


K(d)(w,x) = (akw- Kw) +k (5 x); e «(2 x)): 


假定 入 > 1, di(A) 是 问题 MAx + kOx/(1 4+ 0) = dy 的 特征 值 ， 由 变 分 原理 ， 显 然 d;( 入 ) 关 
于 入 在 [1, 十 oo0) 是 递减 的 ， 且 满足 


—00 = +++ < dg(A) < da(à) < di (à), di (1) = do. 


若 入 是 K(d) 的 特征 值 ， 对 应 的 特征 函数 为 (w, x)， 则 易 得 和 > 1 是 天 (d) 的 特征 值 当 且 仅 当 
存在 某 个 i， 使 得 d = qd;( 入 )。 

定理 3 在 定理 1 的 假设 下 ， 由 (do; 9,0) 处 产生 的 局 部 分 支 可 以 延 拓 成 整体 分 支 ， 且 分 歧 曲 
线 沿 参数 d 到 达 无 穷 远 处 。 

证 明 假设 d > do W 


d > do=di(l) >d,(A), i=1,2,---, 
WH KA 没有 大 于 1 的 特征 值 ， 从 而 idex(T(d,)) = (一 1)? = 1。 假 设 
do = di(1) >d> max{d2(1), 0}, 


Mild > dj(A), i = 2,3,0 由 于 lim di (A) = 一 00， 故 存在 惟一 的 入 > 1 使 得 d = di (Ao), Al 
此 K (d) 存在 大 于 1 的 惟一 特征 值 Xo。 
下 面 说 明 Ao 是 KK(q) 的 单 重 特征 值 ， 只 需 证 明 


N (Aol — K(d))NR(NI- K(d)) = {0} 


即 可 中。 从 而 index(T(d,:)) = (—1)! = -1. 
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由 文献 [6,8] 中 全 局 分 歧 理 论 可 知 ， 在 R+ x 其中 ，G(d;w,x) 存 在 过 点 (do;0,0) 的 连通 分 
支 Co， 且 在 (do;0,0) 附近 G(d;w,x) = 0 的 所 有 正解 都 在 该 分 歧 曲 线 上 。 令 Ci 表示 其 最 大 连通 
分 文 ， 并 定义 


Cı = Co — {d(s), s(y1 + 9(s)), s(%1 + ¥(s)) :—6 < s < 0}, 
则 在 分 歧 点 (do; 0, 0) 附近 C1 包含 曲线 
{d(s), s(el+Toe(s)),s(Wi +W(s)) :0<s< ô}. 


SC = {(d50 —w,x) : (dw, x) ec ， 则 C 是 (3) 从 (9,0) 分 歧 出 来 的 正解 ， 且 在 (do; 0,0) 的 
小 邻 域内 C CP, FH 


P= {(d;u,v) € R* x X:u,v>0,2€Q, =, St < 0, ze 00), 
并 且 分 支 C 满 足下 列 条 件 之 一 。 

(i) C 连 接 分 歧 点 (do;9,0) 到 (6d;9,0)， 其 中 4d 关 do 且 I 一 KK(d) 不 可 道 ; 

(i) C 连 接 (do; 6,0) 沿 着 参数 延伸 至 无 穷 

(iii) C 包 含 形 如 (d;9 — u, v)(d;0 +u, —v) HA, Hu £0, v #0. 

下 证 (ii) IZ. WC 一 {(do;6,0)} Z P, REE AR (d;7,0) € C 一 {(do;9,0)} c OPA 
FRR (dntinsUn) €E COP, f&n 一 00 时 (dn; un,vn) > (3,0. HRA HREM, T= 
ORS = 0. HF (0,0) 4EIBtL, W (G0) = (0,0). Wn 一 colt (dn; un, vn) 一 (d;9,0). 
S Va = Un/|lUnlloos WW Vn 满足 





My EV ee Eo. SER VEO eD. (8) 


—AV, — 
i 1+ Un 1 工 十 nm 


H Le 估计 和 Sobelov 棋 入 定理 ， 必 存在 一 个 收 敏 子 列 ( 仍 记 为 你 ) 使 得 m 一 OM V, > V 
ECM ERZ, HV >0,V 40, ce. Ak, (8) Fen 一 co 得 


k0 > 
AV rgo EV 0, zEQ, V=0, TEDN 


由 极 大 值 原理 得 V > 0。 因 此 = -d,(-k0/(1+0)) = do Gd 关 do 矛盾 。 因 此 C 一 
{(do; 9,0)} C P。 由 引 理 1 中 解 的 先 验 估计 知 ， 存 在 常数 M > 0， 使 得 ||wllw, lulo <M. K 
此 C 在 正 锥 已 内 只 能 沿 参 数 d 趋 于 无 穷 远 处 。 
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Stability of Bifurcation Solutions to an Elliptic Equation for a Class 
of Predator-prey Models 
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Abstract: Under the Dirichlet boundary conditions, the steady-state problem for a predator-prey 
model with predator cannibalism is considered in this paper. By using the spectral analysis and the 
bifurcation theory, the local bifurcation emanting from the semitrivial solution is obtained, and the 
structure of positive solutions is given. At the same time, with the perturbation theorem to linear 
operators, we prove that the local coexistence solution is unconditionally stable. Finally, by using the 
global bifurcation theory, the local bifurcation is extended to the global one, and the direction of the 
global branch is also established. Our results show that when the parameters satisfy certain conditions, 
the predator and the prey will coexist in a given region. 
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